February 28, 2012 1:13 



WSPC - Proceedings Trim Size: 9in x 6in ShuTanaka'KinkiSymposium2011'arXiv 



1 



QUANTUM ANNEALING AND QUANTUM FLUCTUATION 
EFFECT IN FRUSTRATED ISING SYSTEMS 

SHU TANAKA 

Department of Chemistry, University of Tokyo, 
7-3-1, Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan 
E-mail: shu-t@chem.s .u-tokyo . ac.jp 

KYO TAMURA 

Institute for Solid State Physics, University of Tokyo, 
5-1-5, Kashiwanoha, Kashiwa-shi, Chiba, 277-8501, Japan 
E-mail: r.tamura@issp.u-tokyo.ac.jp 



Quantum annealing method has been widely attracted attention in statistical 
physics and information science since it is expected to be a powerful method to 
obtain the best solution of optimization problem as well as simulated annealing. 
The quantum annealing method was incubated in quantum statistical physics. 
This is an alternative method of the simulated annealing which is well-adopted 
for many optimization problems. In the simulated annealing, we obtain a solu- 
tion of optimization problem by decreasing temperature (thermal fluctuation) 
gradually. In the quantum annealing, in contrast, we decrease quantum field 
(quantum fluctuation) gradually and obtain a solution. In this paper we review 
how to implement quantum annealing and show some quantum fluctuation ef- 
fects in frustrated Ising spin systems. 
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1. Introduction 

In information science, it has been an important issue to propose a novel 
method or improve known methods for obtaining the best solution of op- 
timization problems. Optimization problem is to find the state in which 
the real- valued cost /gain function has the minimum/maximum value. In 
general, since there are huge number of elements in optimization problems, 
we cannot obtain the best solution by a naive method such as full search in 
limited time and present resources. This is because the number of states in- 
creases exponentially with the number of elements. Optimization problems 
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are often represented by binary language of Os and Is. Here the minimum 
unit of information is called "bit" . This situation is the same as the Ising 
model. ^ In the Ising model, the minimum unit corresponds to "spin" which 
can have the value +1 or —1. Then, in many cases, optimization problems 
can be represented by the Ising model with random interaction J^. The 
Hamiltonian is represented as 



Here af is an assigned integer ±1 (not operator*^) at the z-th site. If an op- 
timization problem is mapped onto the Ising Hamiltonian given by Eq. (1), 
the cost function of the optimization problem is expressed as the internal 
energy of the mapped Ising spin system. Thus, to find the best solution of 
optimization problem corresponds to obtain the ground state in the mapped 
Ising spin system. 

Ising models with random interactions often have a difficulty of slow 
relaxation at low temperature. In general, energy landscape of complicated 
Ising model prevents the system from relaxation.^"* In order to avoid such 
a difficulty, many methods have been proposed in statistical physics. The 
most famous one is simulated annealing which was proposed by Kirkpatrick 
et al.^'^ In the simulated annealing, we decrease temperature (thermal fluc- 
tuation) gradually and obtain a solution of optimization problem. In order 
to find the best solution, we should sample states according to transition 
probability defined in some way. Transition probability from a state E to 
the other state S' is often defined so as to be proportional to e''^^^^'"^^^ 
where EiT,) is the eigenenergy of the eigenstate S and /3 is an inverse tem- 
perature (3 = T~^. Here we set the Boltzmann constant ks to be unity. At 
low temperature, because of large 0, a transition between states does not 
occur with frequency. On the other hand, at high temperature, the probabil- 
ity distribution is almost flat. Then a transition from a certain state to the 
other state often occurs. By using the effect of thermal fluctuation, we are 
easy to obtain the best solution by decreasing temperature gradually. Ac- 
tually, Kirkpatrick et al. demonstrated finding a stable state of the random 
Ising spin systems. After that, the simulated annealing is widely adopted for 
many kinds of optimization problems in physics, chemistry, and informa- 



''Hcrc wc do not consider quantum system. In quantum systems represents the z- 
component of the Pauh matrix of the i-th site. Even if we consider the quantum version 
of the Hamiltonian given by Eq. (1), there is no off-diagonal elements. In this paper we 
discriminate c-number and matrix (operator) without or with the hat 




(1) 
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tion science since the simulated annealing is easy to implement. Moreover 
it has been proved that the simulated annealing can find the best solution 
definitely when we decrease temperature slow enough J Then the simulated 
annealing has been regarded as guaranteed general method for optimization 

problems. 

Since there are many types of optimization problems, there have been 
proposed corresponding methods in order to obtain the best solution of in- 
dividual optimization problem. This is a style of information engineering. 
However it is also important to construct a general method for optimization 
problems such as the simulated annealing. Based on this point of view, the 
quantum annealing has been also developed. ^^^"^ In the quantum anneal- 
ing, instead of decreasing temperature (thermal fluctuation), wc decrease 
quantum field (quantum fluctuation) gradually and obtain a solution of op- 
timization problem. The quantum annealing method has been expected as 
a powerful tool to obtain the best solution of optimization problem. How- 
ever there are few examples that the efiiciency of quantum annealing is 
worse than that of the simulated annealing. In order to know when to use 
the quantum annealing, it is a significant issue to investigate microscopic 
nature of quantum field response in simple models. 

The organization of the rest of the paper is as follows. In Section 2, 
we review implementation methods of quantum annealing. In Section 3, we 
consider quantum field response of frustrated Ising spin systems. In Section 
4, we conclude this paper briefiy. 

2. Implementation Methods of Quantum Annealing 

In this section, we review concept of quantum annealing and implemen- 
tation methods. The quantum annealing is expected to be an alternative 
method of simulated annealing and a general framework for optimization 
problems as mentioned above. Here we assume that our target optimization 
problem can be expressed by the Ising model given as 

H. = -YlJi,ata^^, (2) 
where af denotes the a-component of the Pauli matrix defined by 

^-iX^' --Co')' 

Equation (2) is a quantum version of the Hamiltonian given by Eq. (1) and a 
diagonal matrix. Our purpose is to find the ground state of the Hamiltonian 
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given by Eq. (2). In order to perform quantum annealing, we introduce an 
additional time-dependent Hamiltonian 'Hq{t) which represents quantum 
fluctuation. Then the total Hamiltonian is expressed as 

H{t) = n, + n^{t). (4) 

Hereafter we refer to "He and Hq(i) as classical Hamiltonian and quantum 
Hamiltonian, respectively. The most typical form of the quantum Hamilto- 
nian for the Ising spin system is transverse field given as 

Hq(i)^-r(t)^af. (5) 

i 

Although there is an ambiguity of 'Hq(t), we concentrate on the case that 
Hq(t) is the form given by Eq. (5) in this paper. 

In the quantum annealing, we gradually decrease quantum field T and 
obtain the final state at F = 0, whereas we gradually decrease temperature 
T and obtain the final state at T = in the simulated annealing (Fig. 1). 
Recently hybrid-type quantum simulated annealing in which we decrease 
temperature and quantum field simultaneously has developed. ^^'^"^ Quan- 
tum annealing can be categorized by three types summarized in Table 1. 
There are some implementation methods of quantum annealing theoreti- 
cally and experimentally. We can divide into two types of implementation 
theoretical methods according to how to treat time-development. The first 
one is a stochastic method which is realized by the Monte Carlo method. 
The Monte Carlo methods have been widely adopted for analysis of equilib- 
rium state of many-body strongly correlated systems in statistical physics. 

Quantum field (r) 



Simulated annealing 

^^^^^^^^^^^^^H— ^ Temperature (T) 

Fig. 1. Schematic diagram of simulated annealing and quantum annealing. The aim 
of annealing methods is to obtain the best solution at the circle (T = F = 0). In 
the simulated annealing, we decrease temperature T gradually from high temperature 
to (nearly) zero temperature. In the quantum annealing, we decrease quantum field F 
gradually from large value to (nearly) zero. 
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Table 1. Quantum annealing is categorized by three types depending 
on how to treat time-development. 



Theoretical methods 


Experiments 


Stochastic methods Deterministic methods 




Quantum Monte Carlo Real-time dynamics 


Artificial lattices 


(Exact diagonalization, 


(Optical lattice) 


timc-dcpcndent DMRG) 




Mean-field type method 





Since in some cases we face on the difficulty of obtaining the equiUbrium 

state by naive Monte Carlo methods, many efficient algorithms have been 
developed. ^^'^^"'^^ Then if we adopt these engineered methods depending 
on the situation, we can use the Monte Carlo method as a powerful tool for 
quantum annealing. 

The other type of theoretical methods is a deterministic methods. We 
can calculate real-time dynamics by solving the time-dependent Schrodinger 
equation directly^"^^'^'"'^^ or performing the time-dependent density matrix 
renormalization group. There is another type of deterministic method, 
which is mean-field type calculation.'^^ Moreover, variational Bayes infer- 
ence which is a useful method in information science can be regarded as a 
method based on mean-field calculation. This method can also treat large- 
scale problems as well as the Monte Carlo method. Originally, the varia- 
tional Bayes inference has been adopted for many optimization problems. 
Recently a quantum annealing version of variational Bayes inference was 
developed. 

Quantum annealing can realize not only theoretically in classical com- 
puter but also experimentally. Artificial lattices such as optical lattice have 
been made in order to simulate strongly correlated fermionic systems and 

spin systems. Since in artificial lattice systems, we can control parame- 
ters in the Hamiltonian, it is expected that these systems become a "quan- 
tum computer" which can find the ground state of complicated systems. 

In this paper we focus on two theoretical implementation methods. In 
Sec. 2.1, we review the quantum annealing method for the transverse Ising 
model by using the quantum Monte Carlo method. In Sec. 2.2, we explain 
the real-time dynamics using the Schrodinger equation for the transverse 
Ising model. 
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2.1. Quantum Monte Carlo method 

In this subsection, we explain how to implement quantum annealing by the 
quantum Monte Carlo method. Here we take the transverse Ising model as 
an example. The quantum Monte Carlo method is a useful tool to obtain 
thermodynamic properties of strongly correlated many-body systems such 
as quantum spin systems, bosonic systems, and fermionic systems. Before 
we review the quantum Monte Carlo method for the transverse Ising model, 
we show implementation method of classical Monte Carlo method for Ising 
spin systems without transverse field given by Eq. (1). The main procedure 
of the classical Monte Carlo method is as follows. 

Step 1 We prepare an initial state''. 

Step 2 We select a spin randomly. 

Step 3 We flip the selected spin with probability defined by some way. 

Step 4 We repeat Step 2 and Step 3 until physical quantities converge. 

Performance of Monte Carlo method depends on how to choice the transi- 
tion probability in Step 3. The simplest definition of the transition proba- 
bility is the Metropolis method. Let S and S' be the state before and after 
spin flip, respectively. In the Metropolis method, we change the state from 
E to E' according to the following probability P(E'|E): 



where -E(E) is the cigenenergy of the eigenstate E. In this rule, when the 
eigenenergy of the candidate state E' is lower than that of the present state, 
the state changes from |E) to |E') definitely. Here |E) represents the direct 
product of A''-spin states as follows: 



where \ai) denotes the state at the i-th site. It should be noted that we 
can easily calculate i?(E) since now we consider the Hamiltonian in which 
there is no off-diagonal element. 

Next we review the quantum Monte Carlo method for the transverse 
Ising model given by 




{E{E') > i?(E)) ' 



(6) 



E) = \ai) \a2) • • • O \aN) 



(7) 




(8) 



''Tlicrc is an ambiguity how to prepare the initial state. For example, we can use both 
random spin configuration and completely polarized configuration. 
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In this model for large number of spins, it is difficult to obtain all eigenen- 
ergies and eigenstates in practice because of off-diagonal elements in the 
Hamiltonian. Then we cannot use the Monte Carlo method directly since we 
should know all eigenvalues and corresponding eigenstates. Wc have to con- 
sider an alternative representation which enables us to treat this model. In 
general, partition function of d-dimensional transverse Ising model (quan- 
tum system) is equivalent to that of (rf+l)-dimensional Ising model without 
transverse field (classical system). This correspondence can be derived by 
the path-integral representation (in other words, the Trotter-Suzuki decom- 
position45'46). 

Then we can use the Monte Carlo method for the transverse Ising 
model by considering path-integral representation. For simplicity, in order 
to demonstrate path-integral representation, we consider the ferromagnetic 
Ising chain with transverse field whose Hamiltonian is given by 

A' N 

n = n, + n^, n, = -jJ2<^t^Ui' ^q = -rl]^r, (9) 

i=l i=l 

where the periodic boundary is assumed: ajv+i = ci. At the temperature 
r(= /3~^), the partition function of this system is expressed as 

-/3CHc+«q) 



(10) 



We cannot obtain the partition function by this direct expression since 
(S|e~^(^<=+^''^|E) is not tractable. Then by using integer m, we decompose 
the matrix exponential e"'^'^^""'"^''^ as follows. 

2\ 



exp [-P{Uc + Hq)] = exp 




Z = 13 (S I e-''(«=+«-) 

E 

= E (^^ 



Then we obtain the partition function by using m: 



P U 



X 


(S2 e- 


X 




X 





Si), 



(11) 



(12) 



where represents the direct product of N spin states as well as |E): 
|E/c) = |cr/c,i) (g) \ak,2) (8) • • • 18) |o-/c,Jv) • (13) 
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The index k represents the coordinate along the Trotter axis. Hereafter we 
refer to m as the Trotter number. Since the classical Hamiltonian He is a 
diagonal matrix, the following relation is satisfied: 



IS.) = C7l 



From the above relation and simple calculation, 

iV X N 

m 



(14) 



E'fe)=exp 



i=l 



Efc+1 ) = 







N 
2 


isinh( 




exp 







- log coth 



N 



\ i=l 



are obtained. Then the partition function is expressed by the following way: 
^1 



Z= lim 

m— >-oo 



„ sinh f I 

2 \ m J 



X exp 

{<Tfc,,=±l} 



N 



(3J 



TO 



lim A 

m— ^oo 



EE 

,1=1 k=i 

exp(-/3'Heff), 



k,i'^k,i+i + 2^ogcot]i 



J2 exp(-/37^eff), (15) 

where A is just a coefficient which is irrelevant for thermodynamic proper- 
ties and the eff'ective Hamiltonian TieS is defined as 



N 



^e. = EE 



J 



TO 



'k,i"k,i+l 



^ log coth a^^a^,^^. 



(16) 



i=l fe=l 

This relation means the partition function of one-dimensional transverse 
Ising model is equivalent to that of two-dimensional Ising model without 
transverse field. Since physical quantities can be calculated by the partition 
function in general, we can obtain thermodynamic properties of transverse 
Ising chain by considering that of two-dimensional classical Ising model. 
It should be noted that the abovementioned derivation does not depend 
on spatial dimension. Then the method can be adopted for general Ising 
model with transverse field. Figure 2 displays schematic spin states of the 
transverse Ising chain for small F and large F at low temperature. When the 
transverse field F is small, spins are almost aligned along both real space 
and the Trotter axis. On the other hand, spins are aligned along only real 
space when the transverse field F is large. This is a nature which comes 
from quantum fiuctuation effect. 
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Trotter axis 



Trotter axis 



tttttttl'ttt 



'i>i't>i'>l'>l^>i'il'tl- 
itttttttttt 



Real space 



Real space 



Fig. 2. Schematic spin configurations of the transverse Ising chain given by Eq. (9) for 
small r (left panel) and large F (right panel). 

So far, we review quantura Monte Carlo niethod as a tool to obtain the 
thermodynamic equilibrium properties. This method can be also regarded 
as a realization method of stochastic dynamics. Then we can use quantum 
Monte Carlo method as a method to perform quantum annealing. In the 
quantum annealing, the transverse field decreases against Monte Carlo step. 
The efficiency of quantum annealing by quantum Monte Carlo method has 
been studied, and the quantum annealing method has been succeeded to 
obtain not so bad solution of some optimization problems. 

2.2. Real-Time Dynamics 

In the previous subsection, we explained implementation method of quan- 
tum annealing by quantum Monte Carlo method. In the method, time- 
evolution is treated as a stochastic process. In this subsection, we review 
a method in which time-evolution is considered as a deterministic process. 
There are a couple of methods which realize deterministic time-evolution. 
In this paper, we focus on real-time dynamics which is derived from the 
Schrodinger equation: 



where {ipit)) and 'H(t) denote time-dependent state and Hamiltonian, re- 
spectively, at time t. In this method, once we prepare a state as the initial 
state, the time-development of the state is determined. Then, how to choice 
the initial state is important. In the sense of quantum annealing, it is a usual 
way to prepare the initial state which can be made easily. In the case of 
the transverse Ising model, fully polarized state is a trivial state at large 
transverse field. The state is expressed as 



ih- m)) = n{t) \m) , 



(17) 



(18) 
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where the state |— )•) is defined as 

K) = i=(|t) + |i)). (19) 

The state |— >■) is an eigenstate of tlic x-componcnt of the Pauli matrix . 
If the energy level of the prepared initial state does not cross with other 
levels for all transverse field strength F, we can obtain the ground state in 
the adiabatic limit. However, in practice, since we decrease the transverse 
field with finite speed, transition to excited levels is inevitable. In order to 
consider such a nonadiabatic transition, we study a single spin system with 
longitudinal and transverse fields. The Hamiltonian is given by 

^sir^gie = -hu'- Va- = ( I J "^^ ) • (20) 

Eigenvalues of this Hamiltonian are 

= ±V/i2 + r2. (21) 

Figure 3 denotes the eigenenergies as functions of longitudinal field h. Then 

the energy difference (energy gap) between two eigenstates is 2\//?"+T^. 
The energy gap takes the minimum value 2F at /i = for F. 




h 



Fig. 3. Eigenenergies of the single spin with longitudinal and transverse magnetic fields. 
The dotted lines and the solid curves indicate the eigenenergies for F = and F = 1, 
respectively. 
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Next we consider a single spin problem with time-dependent longitudi- 
nal field h{t) = vt and the fixed transverse field F: 

Wsingie(t) = -vta^ - = ( ;Q . (22) 

Here we set the down state |4-) as the initial state. This state is the ground 
state of the Hamiltonian in the limit of t — oo. When the transverse 
field r is zero, the state \\.) is the eigenstate of the Hamiltonian. Prom the 

viewpoint of quantum annealing, if the symmetry of the prepared initial 
state is different from that of the ground state at the final time, we cannot 
obtain the ground state at all. We next consider the case for finite transverse 
field r. In this case, the; energy level of the ground state does not cross with 
that of the excited state as stated before. Then the state becomes |t) at 
t — ^ 00 in the adiabatic limit (v — >■ 0) . When we sweep longitudinal magnetic 
field with finite speed v, linear combination of eigenstates is obtained even 
at i — >■ oo. 

The left and right panels in Fig. 4 show time-development of magnetiza- 
tion along the 2;-axis for various sweeping speed with fixed transverse field 
r = 0.2 and for various F with fixed sweeping speed v = 0.1, respectively. 
Here the magnetization along the .z-axis is calculated by 

m^{t) = mw\m)- (23) 

As the sweeping speed becomes slow and/or the energy gap becomes large, 
the state approaches the adiabatic limit of the ground state {m^{t = 
+oo) = -1-1). This is called the Landau-Zener-Stiickelberg transition. "''^"'^^ 
The asymptotic behavior of the transition probability at t ^ -|-oo is given 




Fig. 4. Time-development of the magnetization along the z-axis for various sweeping 
speed under fixed T = 0.2 (left panel) and for various V under fixed sweeping speed 
V = 0.1 (right panel). 
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by 

P.zs^oxp(-i^), (24) 

where AE is the energy gap at h — 0; in this case AE — 2r. The Landau- 
Zener-Stiickelberg transition is adopted not only for single spin problem but 
also many-body quantum systems. Actually, a couple of quantum dynamical 
behaviors can be analyzed by the Landau-Zener-Stiickelberg transition since 
in many cases, the energy level structure can be approximated by that of 
single spin system. Then, in order to consider nonadiabatic transition in 
the quantum annealing, the knowledge of the Landau-Zener-Stiickelberg 
transition is useful. 

3. Quantum Field Response of Frustrated Ising Systems 

In this section we consider transverse field response of frustrated Ising spin 
systems. In order to explain the concept of frustration, we first consider the 
Ising model on a triangle cluster. The Hamiltonian of this model is given 

by 

= -J + alal + a^a^) , = ±1. (25) 

When the coupling constant J is positive/negative, this is called ferromag- 
netic/antiferromagnetic interaction. Figure 5 (a) shows the ground states 
in the case of ferromagnetic interaction. In the ground states, all interac- 
tions are energetically favorable. On the other hand, when the interactions 
are antiferromagnetic (J < 0), in the ground state, there is a bond which 




Fig. 5. Ground state spin configuration of tlie Ising model on a triangle cluster. The 
dotted and solid lines represent ferromagnetic and antiferromagnetic interactions, re- 
spectively, (a) Ferromagnetic case. All interactions are energetically favorable states, (b) 
Antiferromagnetic case. The crosses denote energetically unfavorable states. 
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is energetically unfavorable shown in Fig. 5(b). Such a nature is called 
frustration. Since frustration prevents the system from ordering and makes 
peculiar density of states, unconventional order and characteristic dynamic 
behaviors appear. 

Frustration appears in antiferromagnet on triangle-based lattices such 
as triangular lattice, kagome lattice, and pyrochlore lattice. In random Ising 
spin systems, frustration also appears randomly on the lattice. Since many 
optimization problems can be mapped onto random Ising models as men- 
tioned above, it is an important issue to investigate quantum field response 
of frustrated systems comparing with thermal fluctuation effect. 

In this paper we focus on two examples of frustrated Ising spin systems. 
In Sec. 3.1, we review transverse field effect on fully frustrated systems. In 
Sec. 3.2, we study decorated bond systems in which correlation function 
behaves nonmonotonic against temperature and transverse magnetic field. 

3.1. Order by Disorder Effect in Fully Frustrated Systems 

Here we consider fully frustrated systems where all plaquettes are frus- 
trated. There are macroscopically degenerated ground states in the fully 
frustrated systems. ^^"^"^ Typical configurations of ground state of the anti- 
ferromagnetic system on the triangular lattice are shown in Fig. 6. In Fig. 6, 
we adopt the periodic boundary condition. The dotted boxes in these fig- 
ures indicate "free spin" where the internal field from the nearest neighbor 
sites is zero. The number of free spins is a useful way to represent a char- 
acter of each configuration. From left to right in Fig. 6, the number of free 
spins decreases. Here we first consider simulated annealing for such fully 
frustrated systems. From the principle of equal weight, all distinct degen- 
erated ground states can be obtained with the same probability at T = 0-|- 
in the simulated annealing. 

By the way when we adopt quantum annealing for fully frustrated sys- 
tems, whether does the probability distribution of the obtained ground 
states become flat or not? If the probability distribution becomes biased 
distribution, which arc states selected? From the preceding studies, it is 
well-known that the states which have many free spins tend to appear in the 
quantum annealing. ^^'^^'^"'^^ The left panel in Fig. 6 represents the maxi- 
mum free spin state. Actually the probability of this state is the maximum 
in the adiabatic limit. This is because the transverse fleld is represented as 




(26) 



February 28, 2012 1:13 



WSPC - Proceedings Trim Size: 9in x 6in ShuTanaka'KinkiSymposium2011'arXiv 



14 



Fig. 6. Typical ground states of antiferromagnetic Ising spin system on triangular lat- 
tice. The dotted boxes represent free spin where the internal field from the nearest 
neighbor sites are zero. 



and corresponds to sum of the spin-flip operator. Then the states which 
have large number of free spins appear with high probabihty. Schematic 
picture of the probabihty distributions of the ground state obtained by 
the simulated annealing and the quantum annealing is shown in Fig. 7. 
The maximum free spin states are "ordered states" which are mediated by 
quantum fluctuation effect. Then, the nature is called "order by disorder" 
which is a famous feature in fully frustrated systems. ^^"^'^ The energy of the 
states obtained by simulated annealing with slow schedule is the same as 
that of the states obtained by quantum annealing. Although the maximum 
free spin states are selected in this case, the other type of ground states 
can be selected when we adopt the other type of quantum fluctuation. It 
is an interesting problem to investigate such probability distribution when 
we decrease temperature and transverse field with finite speed. ^® 



(a) 



(b) 



Probability 



Probability 



iiiiiiiiiL LiuiOlL 



Number of free spins 



Number of free spins 



Fig. 7. Schematic picture of probability distributions of the ground state for (a) sim- 
ulated annealing and (b) quantum annealing. The number of free spins represents the 
configuration of ground states. 
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3.2. Decorated Bond Systems 

Next we consider quantum field response of decorated bond systems. Figure 
8 shows a structure of decorated bond systems. The circles and triangles 
in Fig. 8 denote system spins and decorated spins, respectively. Before we 
consider thermodynamic properties of lattice system with decorated bond 
shown in Fig. 8 (b), we first consider temperature dependent correlation 
function between system spins of the decorated bond unit shown in Fig. 8 
(a). The Hamiltonian of the unit is given as 

^unit = - JdirfT^^I - (aj + ) Sf , (jf = ±1, = ±1, (27) 

i=l 

where Jdir and J represent direct coupling between system spins and dec- 
orated bond between a system spin and a decorated spin, respectively. 
The correlation function between system spins can be calculated exactly 

as 

(crj^ cr^ tanhXoff, (28) 
Tr{,je~^^ = ^gifoffaja-^ (-29) 

Xeff = ^ l0gCOsh(2/?J) + /?Jdir, (30) 

where A is an irrelevant factor. Here if J^i^ = and J > 0, the correla- 
tion function is always positive and monotonic decreasing function against 
temperature. On the other hand, when J = and J^,^ < 0, the correlation 
function is always negative and monotonic increasing function against tem- 
perature. Then, by tuning the ratio Jdir and J with keeping J^ir < and 
J > 0, the correlation function (uf erf) behaves non-monotonic as a function 




Fig. 8. (a) Unit of the decorated bond system, (b) Square lattice with decorated bond 
for = 1. 
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of temperature. ^^'^^'^^^''^ In this paper we adopt Jdir = — ^ J and J as an 
energy unit. The temperature dependency of correlation function in this 
case is shown in Fig. 9. It is noted that the signs of correlation function 
and the effective coupling K^s have both plus and minus values depending 

on temperature. 

In order to investigate the reason for such a non-monotonic behavior, 
we consider probability distributions of ferromagnetically correlated states 
(afcrf = +1) and antiferromagnetically correlated states (crfcf = —1). The 
number of ferromagnetically correlated states and that of antiferromag- 
netically correlated states are 2^''+^. Let A^s be the number of decorated 
spins in which directions of these spins are the same as af. The probability 
distributions of ferromagnetically correlated states Pp(-^s) and antiferro- 
magnetically correlated states Paf(-^s) are given by 



PAF(iVs) 



2 cosh (/3^J) 



Nd g-2/3(Wd-2W,)J 



cosh^" {2/3 J) ' 



2 cosh J) \ KJ \2 



(31) 
(32) 




Fig. 9. Correlation function between system spins as a function of temperature. 
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It should be noted that the following relation is obviously satisfied: 

iVd 

^ [Pf{N,) + Paf{N,)] = 1. (33) 

At zero temperature all system spins and decorated spins arc the same 
value. Then the probability distribution is as follows: PpiNa) — 1 and 
otherwise are zero. At the temperature where J2n Pf{^s) = J2n ^af(A^s) 
is satisfied, the correlation function (erf cf) becomes zero. This is the reason 
why the correlation function and effective coupling behave non-monotonic. 

Suppose we consider square lattice system with decorated bond for large 
enough N^. In this case there is temperature region where the absolute 
value of effective coupling exceeds the critical value of the Ising spin sys- 
tem on square lattice Kc — ^ log (l + -\/2) Then, paramagnetic phase 
— >■ antiferromagnetic phase — > paramagnetic phase — > ferromagnetic phase 
appear as temperature decreases and at each phase boundary, phase tran- 
sition takes place. Such successive phase transitions arc called reentrant 
phase transition. Reentrant phase transitions and non-monotonic behavior 
of correlation function are typical nature of frustrated systems. 

So far we showed thermal fluctuation response of decorated Ising spin 
systems. Next we consider transverse field response of this system at zero 
temperature. Then the Hamiltonian is given by 

n = - Jdiraf a| -jJ2{^t+ ^1) - r + 5f + 5^ an , (34) 

i=l \ i=l / 

where sf denotes the a-element of the Pauli matrix of the i-th decorated 
spin. As the previous example, we adopt Jair = — ^J. We calculate the 
correlation function between system spins along the z-axis at zero temper- 
ature as a function of transverse field as shown in Fig. 10. As well as the 
thermal fluctuation, the correlation function behaves non-monotonic as a 
function of transverse fleld. This result indicates there is a similar point 
between thermal fluctuation and quantum fluctuation originated by trans- 
verse field. Suppose we consider the system where reentrant phase transition 
occurs when we decrease temperature. If we adopt the simulated anneal- 
ing method, we face on the difficulty which comes from the critical slowing 
down and the simulated annealing is not useful. Whenever we adopt the 
quantum annealing using time-dependent transverse field, the system also 
exhibits reentrant phase transition. Then we should consider other type of 
quantum fluctuation effect in order to avoid the occurrence of the phase 
transition. 
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Fig. 10. Correlation function between the system spins as a function of transverse field. 
4. Conclusion and Future Perspective 

In this paper, wc reviewed how to implement quantum anneahng foeus- 
ing on the methods based on quantum Monte Carlo method and time- 
development Schrodinger equation. The quantum annealing method is ex- 
pected to be a powerful tool to obtain the best solution of optimization 
problems. This method is a general method and can be implemented easily. 
In this method we use quantum fluctuation to find the stable state. In order 
to put into practical use, it is important to investigate quantum fluctuation 
effects for optimization problems. Then we considered quantum fluctuation 
effect for frustrated systems by taking the transverse Ising model as an 
example, since the Ising model which represents optimization problem is 
a frustrated system in general. We explained differences and similarities 
between thermal fluctuation and quantum fluctuation. By building up the 
organized knowledge of the thermal fluctuation and the quantum fluctua- 
tion, the quantum annealing method develops into a truly useful algorithm. 
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